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RESUMEN 

Los Estallidos de Rayos Gamma (ERG's) largos son las explosiones electromagneticas mas potentes del Uni- 
verso, asociadas a la muerte de estrellas masivas. Como tales, los ERG's son trazadores potenciales de la 
evolution de la formation estelar cosmica, la metalicidad y la funcion initial de masa. Los ERG's tambicn 
han probado ser atractivos como indicadores de distancia cosmologicos, lo cual abre una oportunidad unica de 
constrenir la historia de expansion cosmica hasta z ~ 5 — 6. Se presenta aqui una resena sobre ambos temas. 

ABSTRACT 

Long Gamma-Ray Bursts (GRBs) are the brightest electromagnetic explosions in the Universe, associated to 
the death of massive stars. As such, GRBs are potential tracers of the evolution of the cosmic massive star 
formation, metallicity, and Initial Mass Function. GRBs also proved to be appealing cosmological distance 
indicators. This opens a unique opportunity to constrain the cosmic expansion history up to redshifts 5-6. A 
brief review on both subjects is presented here. 
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1. INTRODUCTION 

Detected as brief (~ 0.01 — 100 s), intense flashes 
of 7-rays (mostly sub-MeV), Gamma-ray Busts 
(GRBs) are the brightest electromagnetic explosions 
in the Universe. The power emitted by GRBs in 
electromagnetic form can reach luminosities up to 
L ~ 10 52 - 10 53 erg s -1 , while AGNs can have 
L ~ 10 48 erg s _1 (but for long times), and Su- 
pernovae can have L ~ 10 45 erg s _1 for the first 
hundreds of seconds after the explosion. The short 
variability timescales of the 7— ray emission, suggest 
already very small dimensions for the sources, of the 
order of tens of kilometers, typical of stellar black 
holes or neutron stars. Several pieces of evidence 
indeed show that GRBs arc associated with cata- 
clysmic stellar events, and that the 7— ray emission 
comes from highly relativistic collimatcd outflows. 
The typical bulk Lorentz factor for the jets is T ~ 
300. Thus, GRBs are true cosmic laboratories for 
the study of relativistic, magneto-hydrodynamical, 
and high energy processes (for recent reviews on the 
GRB physics see e.g., Zhang & Meszaros 2004; Piran 
2005; Meszaros 2006). 

Furthermore, GRBs and their afterglows are of 
great interest for studies related to stellar astro- 
physics, the interstellar and intcrgalactic medium, 
and most important, they reveal themselves as 
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unique probes of the high rcdshift Universe. In the 
last 3 years, on average ~ 1 paper per day is pub- 
lished on GRBs in refcreed journals and ~ 4 per day 
in non-refereed publications. 

GRBs are divided into two main groups which, 
following the notation of Zhang (2007), we will call 
Type I and Type II. The former (« 1/3) have short 
7— ray durations (< 2 s) and hard spectrum; it is 
conjectured that they result from binary mergers of 
compact stellar objects (NS-NS or NS-BH). The lat- 
ter (« 2/3) have durations larger than 2 s and their 
7— ray spectra tend to be softer. The observations 
show (e.g., Hjorth ct al. 2003; Stanek et al. 2003) 
that these GRBs result from the collapse of rapidly 
rotating cores of low-metallicity stars more massive 
than about 25 M ('collapsar' scenario). 

A breakthrough in the GRB field happened a 
decade ago: thanks to their relatively long duration, 
for some Type II GRBs it was possible to resolve 
and detect the afterglow at softer (X-ray, optical, IR 
and radio) energies. This allowed to measure spec- 
tral lines and/or to identify the host galaxy; hence 
the redshifts, z, could be determined. Up to Decem- 
ber 2007, there were around 100 GRBs with secure 
z measurements (in ~ 90% of the cases by the after- 
glow and ~ 10% of the cases by the host galaxy). 
More than 60% of the z determinations were ob- 
tained during the last 3 years with the dedicated 
Swift satellite (Gehrels et al. 2004), which allowed 
also to discover the afterglows for some Type I events 
(7 up to date). No doubt, GRBs are the brightest 
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transient cosmological events measurable. In Fig. 1 
the z distribution of Swift GRBs is shown: it ranges 
from z = 0.033 to 6.29, with an average z of 2.2. In- 
direct estimates suggest, however, that many of the 
observed Type II GRBs could be produced at z's 
larger than 6 (e.g., Bromm & Loeb 2005). 

Summarizing: type II GRBs are extremely pow- 
erful explosions associated to the collapse of short- 
lived massive stars, with peak emission at sub-MeV 
energies, where dust extinction is not an issue. Be- 
sides, the 7— ray spectra of type II GRBs are such 
that the K correction is small or even negative. 
Thus, the fluxes of these events can be detected even- 
tually from any z. All these properties convert type 
II GRBs in cosmological beacons, which can help us 
look back in time. We will review here advances 
along this line in two directions: GRBfjf] as tracers 
of the history of global massive star formation rate 
(SFR; §2), and GRBs as cosmic rulers able to help 
us in constraining the expansion history of the Uni- 
verse (§3). In §4, perspectives for future work will 
be discussed. 

2. GRBS AS TRACERS OF THE GLOBAL 
MASSIVE STAR FORMATION RATE 

The death rate of massive short-lived stars re- 
sembles their formation rate. Thus, the GRB for- 
mation rate (GFR) can be used as a potential tracer 
of the massive SFR in the Universe. Current obser- 
vations allow to construct the (yet very incomplete) 
GRB z distribution per unit of time, N(z). This ob- 
servable distribution is connected to the history of 
the intrinsic GFR (per unit of comoving volume), 
^GRb(z), through: 



L-V(Psei§o"zj 



N(z) 



dV/dz F(z) . 
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where dV/dz is the comoving volume element, (1 + 
z)^ 1 accounts for the time dilation, f b s stands for 
the detector exposure factor and the average GRB 
beaming, and F(z) = 4'fiux{z)4>z(z) takes into ac- 
count several selection effects. F(z) can be under- 
stood as the probability to detect the burst and its 
afterglow, and to measure its z from the afterglow; 
it comprises two classes of effects: 

• the flux-limited selection function, 



<t>flux(z) 



L iso (P m i„,z) 



/(Liso)dLj SO , 



(2) 



which depends on the detector flux threshold 
P m i n and on the GRB luminosity function (LF), 




'Hereafter GRB refers to a type II burst. 



Fig. 1. Normalized distribution of 72 spectroscopically 
measured z's obtained with Swift (solid line) and 220 
pseudo-redshifts obtained with the Li so -variability cor- 
relation (dashed line). 

/(Li so )dLi SO (Lj SO is the equivalent isotropic beam 
uncorrected- GRB bolomctric luminosity). 

• the selection function related to the observabil- 
ity of the afterglow and its z determination, (f> z (z). 
This function is very uncertain and comprises spec- 
troscopic and photometric selection effects (e.g., ob- 
servability of lines in the spectral range and in the 
presence of night-sky lines, afterglow flux-limited se- 
lection, different systematics in the emission or ab- 
sorption line technique used to determine z), as well 
as potential astrophysical effects (e.g., obscuration, 
evolving dust extinction, a bias in the GRB host 
galaxy population). An evidence that these selec- 
tion effects are severe is the fact that only ~ 1/3 
of the Swift GRBs have reliable z determinations 
in spite of the quick afterglow localization and the 
great effort in follow-up observations. 

2.1. Direct inferences, yet limited 

A direct joint determination of yicRB^) and 
/(Li so )dLi SO would be possible from a populated 
Liso — z diagram. However, the situation is compli- 
cated mainly by some of the uncertain component 
biases present in the (j> z (z) function (e.g., Bloom 
2003; Gou et al. 2004; Fiore et al. 2007; Coward 
et al. 2007). In view of this and the yet low popu- 
lation statistics, it is not feasible to attempt a clean 
reconstruction of ugrb {z) from current data. Re- 
cent works that make use of the Swift data start 
from some assumptions and limit themselves to ex- 
plore only statistical compatibilities of the observed 
N(z) distribution with models of nGRB(^) assumed 
proportional to the global SFR history inferred from 
extragalactic studies, Psfr(z). A general conclu- 
sion is that N(z) implies an enhanced tt-grb^) at 
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high redshifts with respect to psfr(z) (Kistler et al. 
2008; see also Le & Dermer 2007; Guetta & Piran 
2007). It should be noted that some pieces of evi- 
dence suggest that the Swift sample of GRBs with z 
determined is a fair sample of the real high— z GRB 
population (Fiore et al. 2007). 

On the other hand, Coward et al. (2007) con- 
cluded that the observed Swift z— distribution, 
where the number increasing from z rs to z ~ 1.5 
is modest, implies a bias in the afterglow observabil- 
ity such that at z = — 1 it works inversely propor- 
tional to the global SFR history. It could be that the 
Swift sample with z is biased against low-z GRBs, 
probably due to the enhanced extinction associated 
with the prolific SFR at z ~ 1 (Fiore et al. 2007; 
Coward et al. 2007). 

2.2. Indirect inferences, encouraging results 

Alternatively to the observed L; so — z diagram, 
there are other methods for inferring the GFR his- 
tory from observations but in a model-dependent 
way. We will mention three methods: 

(1) The most extensive GRB observational 
database is the CGRO-B ATSE peak flux P distri- 
bution, N(P), for ~ 3000 bursts. The N(P) distri- 
bution (corrected by the exposure factor) is the re- 
sult of 4 physical ingredients: ?t.grb(z), /(Li so )dL iso , 
the jet opening angle distribution, f(6j), and the 
volumetric factor given by the cosmological model. 
Therefore, the inference of jigrb {z) from the ob- 
served N(P) is a highly degenerated problem. The 
adequate introduction of complementary observa- 
tional information helps to overcome partially the 
degeneracies; for example, the £" p k distribution 
(Daigne et al. 2006) or the 6j distribution (Le & 
Dermer 2007). 

(2) The Ljso — z diagram can be (indirectly) ob- 
tained for large data-sets of GRBs without z mea- 
sured by using empirical correlations that involve 
Lj S0 . For example, Lloyd-Ronning et al. (2002) 
and Yonetoku et al. (2004) inferred ?t.grb(z) by 
applying the Lj so -variability correlation (Fenimorc 
& Ramirez-Ruiz 2000) to 220 BATSE GRBs, and 
the L iso - Sparest relation to 689 BATSE GRBs, re- 
spectively. This method relies on the certainty and 
accuracy of the used empirical correlation. 

(3) Given an adequate parametrization for 
^GRbO 2 ) an d /(Liso)dLi so , their free parameters can 
be efficiently constrained by fitting models jointly 
to both the observed N(P) distribution and the 
N(z) distribution inferred as in (2) or as in §§2.1 af- 
ter selection effects correction. So far, this method is 
the most powerful. A description of and the results 



from this method (Firmani et al. 2004; hereafter 
FAGT) are as follows. 

The method. The observed N(P) and 
N(z) distributions are modeled by seeding at each z 
a large number of GRBs with a given rate, iigrb(z), 
and LF, /(Li so )dL; so , and then by propagating the 
flux of each source to z = 0. L; so in the rest frame is 
defined as L iso = f^^y eV ES{E)dE, where S(E) is 
the Band (Band et al. 1993) energy spectrum. The 
break energy at rest, Eb, is assumed either constant 
(512 keV) or dependent on L iso according to the 
"Yonetoku relation" (Yonetoku et al. 2004). FAGT 
explored two models for /(Li so )dL; so , the single and 
double power laws (SPL and DPL, respectively), and 
two cases, one where /(Li so )dLi so is constant in time, 
and another one where Lj so evolves as (1 + z) s . The 
function ?tgrb(z) was modeled as 

^grb(z) =K x e S FR(z;a,b) x T) z>2 (z;c), (3) 

where esFR^; a, b) is a bi-parametric function pro- 
posed by Porciani & Madau (2000) for fitting the 
observed SFR history, and r] z> 2(z;c) allows to con- 
trol the growth or decline of hgrb (z) at z > 2. 
The idea is to constrain the parameters of the LF (2 
and 3 for the SPL and DPL, and 1 more if evolu- 
tion is allowed ) and of hgrb {%) (4 parameters) by 
applying a joint fit of the model predictions to the 
N(P) and N(z) distributions. The N(P) distri- 
bution for more than 3000 GRBs compiled by Stern 
et al. (2002) and the z distribution inferred for 220 
GRBs in Lloyd-Ronning et al. (2002; see above) 
were used. 

Results. FAGT obtained that an evolving LF is 
preferred in all the analyzed cases (with an optimal 
value of 5 = 1 ± 0.2). For non evolving LFs, the 
predicted N(P) distributions have systematically 
an excess of GRBs at the bright end, and the peak 
of the N(z) distribution is shifted to higher z's (see 
their Fig. 1). Only a weak preference was found 
for the SPL model over the DPL one. On the other 
hand, the best fits not only imply an evolving LF 
but also GFRs histories that steeply increase (by a 
factor of ~ 30) from z = to z w 2 and then continue 
increasing gently up to z ~ 10 as (1 + z) 14 and (1+z) 
for the SPL and DPL LFs, respectively. 

Figure 2 reproduces from FAGT the SFR history 
as traced by the GFR history (dot-shaded region) 
under an opportune normalization and assuming a 
constant initial mass function (IMF). Both the best 
SPL and DPL LF models with evolution (S = 1) 
and their la uncertainties are taken into account 
by the region. Symbols are the SFR traced by the 
rest-frame UV luminosity and corrected for dust ob- 
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Fig. 2. The cosmic SFR histories as traced by GRBs 
(dotted region) and by rest-UV luminosity (symbols and 
solid line, the latter from the recent fit to several obser- 
vations by Hopkins & Beacom 2006); see text. 

scuration (Giavalisco et al. 2004), while the solid 
line is a more recent piecewise fit to a large compi- 
lation and dust correction of SFR observations by 
Hopkins & Beacom (2006). The main difference be- 
tween the SFR histories inferred from GRBs and 
from extragalactic observations is the enhancement 
(that increases with z) of the former with respect 
to the latter after z = 1 — 2. Interesting enough, 
the Swift based studies mentioned in §§2.1, while in 
much less detail, attain a similar conclusion. 

Implications. Why the GFR evolution might 
be enhanced with respect to the rest-UV traced SFR 
evolution at high redshifts? Let us discuss the pros 
and cons of some of the possibilities. 

• Since the 7—rays do not suffer absorption, even 
in very dense molecular regions, GRBs are expected 
to trace the massive SFR of any galaxy/region in 
the Universe, for example of the dust-enshrouded ac- 
tively star-forming galaxies at high redshifts. How- 
ever, there are pieces of evidence against GRBs being 
located in this kind of galaxies (e.g., Le Floc'h et al. 
2007). 

• Most of the GRB host galaxies are faint, low- 
mass, low-metallicity star-forming galaxies (e.g., 
Fynbo et al. 2003; Stanek et al. 2006) and the GRB- 
to-SN ratio might increase significantly at lower 
metallicities (Yoon et al. 2006). Therefore, the GFR 
could be tracing the SFR of a biased population of 
low-metallicity galaxies. Nuza et al. (2007), by 
means of cosmological simulations, have found in- 
deed that the GRB host properties are reproduced 
if GRBs form from low mctallicty progenitor stars. 
However, recent observational studies are showing 



that the neutral ISM around GRBs is not metal poor 
and is enriched by dust (see Savaglio 2007). Fur- 
thermore, GRB hosts should not to be special, but 
normal, faint, star-forming galaxies (the most abun- 
dant), detected at any z just because a GRB event 
has occurred. 

• If GRBs are produced in binary systems and 
the probability of interloper-catalyzed binary merg- 
ers in dense star clusters (where most GRBs appear 
to occur, Fruchter et al. 2006) increases with z, then 
an enhancement in the GFR is expected (Kistlcr 
et al. 2007). This scenario requires quantitative 
calculations. Interesting enough, FAGT, and more 
recently Bogomazov, Lipunov & Tutukov (2007), 
have calculated the Galaxy production rates of Wolf- 
Rayet stars in very close binaries (progenitors of 
rapidly rotating -Kerr- black holes), and found an 
approximate agreement with the estimated rates of 
type II GRBs in galaxies. 

• An IMF that becomes increasingly top-heavy 
at higher z, would increase the relative number of 
massive stars produced. Since the rest-UV luminos- 
ity already traces massive SFR, the IMF should be 
biased to very massive stars (GRB progenitors) so 
that an evolving IMF could explain the difference 
seen in Fig. 2. 

The above discussion is an example of the large 
potential of GRB studies for exploring such key is- 
sues in astronomy as SF at high z, and the IMF and 
chemical evolution in galaxies. 

Regarding low z's, Fig. 2 shows that the GFR de- 
cays with time closely proportional to the SFR since 
z ss 1. This feature does not seem to be observed in 
the Swift z distribution, which would imply a dom- 
inant rcdshfit bias, 4> z (z) as suggested by Coward et 
al. (2007). In the rcdshift range from z = 1 to z = 2, 
the SFR attains a broad maximum, while the GFR 
keeps increasing rapidly. The Swift z distribution 
would suggest the opposite. In fact, as shown in Fig. 
1, the Swift and the inferred from the L— variability 
relation N(z) distributions differ quite substantially 
up tozw2 (rather than the absolute values, the im- 
portant comparison between the two histograms is 
at the level of the fractional changes with z) . As re- 
marked by Fiore et al. (2007), the present samples 
with z determinations seem to be largely incomplete, 
especially at z ^ 2. 

3. GRBS AS TRACERS OF THE UNIVERSE 
EXPANSION RATE HISTORY 

The energetics of GRBs spans 3-4 orders of 
magnitude; at first sight GRBs are all but stan- 
dard candles. A breakthrough in the field hap- 
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pened after the discovery of a tight correlation be- 
tween the collimation corrected bolometric energy, 
E-y = E- lso (l — cos9j) and the prompt peak energy 
in the vf v spectrum, E p ^ (Ghirlanda et al. 2004a). 
This and any other similar correlations allow to stan- 
dardize the GRB energetics for using GRBs as dis- 
tance indicators in the Hubble diagram (HD). Such 
an endeavor, however, is not trivial. 

The first conceptual problem is that most of the 
GRBs with z measured are at cosmological distances 
(see Fig. 1). Therefore, the given correlation can not 
be calibrated locally; to establish the correlation, a 
cosmology should be assumed, but the cosmological 
parameters are just what we pretend to constrain 
in the HD. By using statistical approaches, this cir- 
cularity problem can be treated in order to get op- 
timal constraints for the explored cosmological pa- 
rameters. The idea is to use the best-fitted corre- 
lation (smallest scatter) for each cosmology and to 
find which cosmology produces the smallest \ 2 m the 
HD, constructed by applying the corresponding cor- 
relation. A powerful Bayessian-like method to carry 
out this undertaking has been introduced by Firmani 
et al. (2005,2007). Other groups have developed al- 
ternative variants (see e.g., Xu, Dai & Liang 2005; 
Schaefer 2007; Li et al. 2007). On the other hand, 
Ghirlanda et al. (2006) have shown that in order 
to calibrate for example the "Ghirlanda" correlation 
it is enough to have a dozen of GRBs in a narrow 
(Az/z ~ 0.1) rcdshift bin, something that will be 
possible in the near future. 

Results. The first cosmological results obtained 
by using the "Ghirlanda" correlation were encourag- 
ing: they provided a test of the accelerated expan- 
sion independent from the SN la studies (Ghirlanda 
et al. 2004b; Firmani et al. 2005; see also Dai, 
Liang & Xu 2004). The addition of GRBs to the 
SN probe, reduces the confidence levels of the con- 
strained cosmological parameters. A drawback of the 
"Ghirlanda" correlation is that it needs to establish 
expensive follow-up observations: 8j is determined 
from the achromatic break time, ibreak, hi the after- 
glow light curve. 

Firmani et al. (2006a) discovered a tight correla- 
tion among three prompt emission GRB parameters: 
Li SO , -Bpk, and duration t45. Some of the cosmological 
constraints obtained by using the "Firmani" corre- 
lation are plotted in Fig. 3. We show in these plots 
how the confidence levels given by the SNLS survey 
(Astier ct al. 2006) are improved when the GRB con- 
straints are added. Right panel of Fig. 3 shows the 
HD for 117 SNLS supernovae and 19 GRBs by us- 
ing the 'vanilla' ACDM cosmology (solid line) , which 



provides a good fit to the observations. GRBs are a 
natural extension of SNIa to high z'. In the bot- 
tom panel the residuals to the assumed cosmology 
are plotted; the averages and its uncertainties are 
0.15 ±0.01 mag and 0.26 ±0.05 mag for the SNe and 
GRBs, respectively 

Our results in general (constraining only two pa- 
rameters at the same time) showed that the flat 
ACDM cosmology is consistent at the la level with 
the HD of GRBs (and GRBs+SNe) up to z = 4.5. 

The cosmography with GRBs opens a valuable 
window for exploring the expansion history of the 
Universe to z > 2 (more than ~10 Gyr ago), where 
SNIa are practically impossible to observe. Besides, 
GRBs offer some important advantages for cosmo- 
graphic studies. (1) GRBs are not affected by dust 
extinction. (2) The luminous distance dr, in the HD 
is a cumulative quantity with z, so that the differ- 
ences among different cosmological models become 
larger at higher z's. Thus, data-points in the HD at 
high z's highly discriminate the models, even if the 
uncertainties are large. (3) Each point in the HD 
translates into a different curve in the space of the 
cosmological parameters (degeneration). The wider 
in z is the data sample, the less elongated along one 
curve (less degeneration) will be the parameter con- 
fidence levels (Firmani et al. 2007). 

Caveats. GRB cosmography is in its infancy 
and of course there should be several caveats as was 
discussed in the literature (e.g., Friedman & Bloom 
2005). For example, it was argued that when using 
the "Ghirlanda" relation, the results are strongly de- 
pendent on the assumption about the density distri- 
bution of the circumbust medium (the dependence of 
6j on tbroak changes slightly with the distribution as- 
sumed, and the parameters of one or another distri- 
bution are included in the calculation of 6j ) . Despite 
that the "Ghirlanda" correlations are different in one 
case or another, from the point of view of cosmogra- 
phy, the results are very similar (Nava et al. 2006). 
Furthermore, it was shown that an empirical corre- 
lation among Ej SO , E p ^, and ibrcak holds (avoiding 
then the assumption of the circumbust density dis- 
tribution), which gives cosmographic results similar 
to those obtained with the "Ghirlanda" correlations 
(Liang & Zhang 2006). Recent Swift observations 
have shown that the A— ray afterglow light curve is 
more complex than previously though and its break 
time tends to be different from the one inferred in 
the optical bands. Several pieces of evidence suggest 
that the A— ray and optical components come from 
different emitting regions; therefore, the requirement 
that the optical tbrcak should be compatible with the 
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Fig. 3. Left: Contours at 68.3% CL's on the (f2 m ,OA) plane from the SNIa HD using our Bayesian approach to circumvent 
the circularity problem (solid green line) and from the combined SNIa+GRB HD (red shaded region) . Middle: Same as 
left panel but on the (fl m , wo) plane assuming flatness; wo is the equation of state index, assumed to be constant. Rigth: 
117 SNIa (red) and 19 GRB (blue) data-points in the HD for the flat ACDM model. The residuals of the data-points 
minus the ACDM model are shown in the bottom panel. 



X— ray one should be relaxed (Nava et al. 2007). 

Concerning the "Firmani" correlation, it was 
established for prompt 7— ray emission quantities 
alone. Therefore, it is model independent and does 
not requires follow-up observations. Some of the 
potential difficulties mentioned in general for the 
methods of standardizing the energetics of GRBs are 
the systematics and outliers in the correlations, the 
gravitational lensing, the possible evolution of GRB 
properties, and the lack of a physical interpretation 
of the correlations. We refer the reader to Ghisellini 
(2007), Firmani et al. (2007), and Ghirlanda (2007) 
for discussions on these caveats. In our opinion, the 
last problem is the most challenging. The possibil- 
ity of evolution is also real (FAGT; Li 2007), but is 
most likely that this happens at the level of the over- 
all population and not in what concerns the internal 
emission mechanisms, which control the spectral en- 
ergy relations. 

Finally, it should be said that the current samples 
of usable GRBs for the correlations involving ener- 
getics are still small. As the samples will increase in 
number, a better treatment of systematics and selec- 
tion effects will be possible. On the other hand, a 
wide spectral 7— ray coverage is necessary in order to 
obtain reliable correlations involving the energetics 
of GRBs. Unfortunately, the Swift BAT detector 
has a too narrow spectral coverage. 

In our view, the use of GRBs as cosmological 
distance indicators has not been sufficiently appreci- 
ated by the astronomical community. We are aware 
of the difficulties of the method, but it should be 
considered that GRBs offer a unique possibility to 



constrain the expansion history of the Universe at 
z > 2, and in a way that simply extends (and comple- 
ments) the method based on SNe type la. Perhaps 
we are now in a similar situation as the SN la workers 
in the early 90 's, when the astronomical community 
used to react skeptically to their proposals. How- 
ever, the effort is worth it, because the cosmography 
with GRBs may offer valuable information to unveil 
the properties of what we call Dark Energy, the big 
mystery of cosmology. This mystery stimulates now 
the frontiers of physics to move in the direction of 
exploring new elements of high energy physics, the 
unification of gravity and quantum physics, gravity 
beyond Einstein relativity, and extra dimensions. 

4. OUTLOOK 

The inference of GFR history and its comparison 
with the cosmic SFR history will highly benefit from 
the growth of the samples of GRBs with accurate z 
determination. However, as discussed in §§2.1, the 
selection effects that plague these samples are a chal- 
lenging issue. Therefore, the indirect methods for in- 
ferring the GFR history (§§2.3) should be improved 
in parallel. These methods would largely benefit 
if large samples of GRBs with pseudo-redshifts are 
constructed. The use of correlations among prompt 
7— ray quantities alone (e.g., the "Firmani" correla- 
tion) is the best way for this aim. 

The improvement of the known GRB tight cor- 
relations or the discovery of new ones, is not only 
important for GFR studies but also for cosmogra- 
phy at high z's as shown in §3. The increasing of the 
observational data is mandatory in this endeavor. 
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With a sample of ~ 150 GRBs, on one hand, the 
empirical correlations might already be calibrated in 
a small z bin, and in the other one, the HD would be- 
come highly populated as to get tight constraints on 
the cosmological parameters. However, it should be 
noted that for the correlations that involve energetics 
and spectral information, a wide spectral coverage is 
necessary, something in which the Swift detector 
fails. The hope is in feature missions as GLAST. 

VA thanks the organizing committee for the in- 
vitation and for a job very well done. Support for 
this work was provided by PAPIIT-UNAM grant 
IN107706-3. 
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